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ABSTRACT 

We report the direct detection of a binary/disk system towards the high- 



> 

■ mass (proto) stellar object IRAS20126+4104 at infrared wavengths. The presence 

of a multiple system had been indicated by the precession of the outflow and 



the double jet system detected earlier at cm-wavelengths. Our new K, L' & 
M' band infrared images obtained with the UKIRT under exceptional seeing 



conditions on Mauna Kea are able to resolve the central source for the first 



time, and we identify two objects separated by ~ 0.5" (850 AU). The K and L' 
£j . images also uncover features characteristic of a nearly edge-on disk, similar to 

many low mass protostars with disks: two emission regions oriented along an 
outflow axis and separated by a dark lane. The peaks of the L' & M' band and 
mm-wavelength emission are on the dark lane, presumably locating the primary 
young star. The thickness of the disk is measured to be ~ 850 AU for radii < 
1000 AU. Approximate limits on the NIR magnitudes of the two young stars 
indicate a high-mass system, although with much uncertainty. These results are 
a demonstration of the high-mass nature of the system, and the similarities of the 
star-formation process in the low-mass and high-mass regimes viz. the presence of 
a disk-accretion stage. The companion is located along the dark lane, consistent 
with it being in the equatorial/disk plane, indicating a disk-accretion setting for 
massive, multiple, star-formation. 

Subject headings: stars: formation - stars: massive - stars: binary - circumstellar 
matter - ISM: individual (IRAS20126+4104) - infrared: ISM 



1 emails: tksridha@cfa.harvard.edu, Stewart.Williams@manchester.ac.uk, Gary.Fuller@manchester.ac.uk 



- 2 - 



1. Introduction 

The study of high-mass star-formation has attracted much recent attention. Molinari et 
al (1996, 2002) and Sridharan et al (2002) have reported systematic studies of candidate high- 
mass proto-stellar objects (HMPOs), leading to a number of follow-ups. A key unresolved 
question is: how similar are the high- and low-mass star-formation processes? Do high- 
mass stars go through a disk-accretion phase, which is a central idea in the current, widely 
accepted scenario for low-mass star-formation? Sufficient observational data is beginning to 
be available for the HMPOs to answer these questions. The outflow phase has been shown 
to be as prevalent towards the HMPOs as towards lower-mass objects, and by implication an 
accretion disk can be inferred (Zhang et al 2001, Beuther et al 2002b). While there is good 
evidence for rotating disks of radius ~ several 1000 AU (Cesaroni et al 1997, 2005, Zhang, 
Hunter & Sridharan, 1998, Beltran et al 2005), outflows presumably require a compact disk 
and direct evidence for their presence is rare (Patel et al., 2005). Compact disks have been 
proposed previously, although not directly imaged (Cesaroni et al 1997 & 1999; Shepherd & 
Kurtz 1999, Shepherd, Claussen & Kurtz, 2001). 

Among HMPO candidates, IRAS20126+4104 is particularly well studied (Cesaroni et 
al 1997, 1999, 2005, Zhang et al 1999, Yao et al, 2000, Shepherd et al 2000). It has a 
luminosity of ~ 10 4 L Q at a distance of 1.7 kpc. Recent observations suggest a mass of ~ 7 
M Q (Cesaroni et al 2005), significantly lower than previous estimates of 20-24 M (Cesaroni 
et al 1999, Zhang, Hunter & Sridharan, 1998). It harbours water, hydroxyl and methanol 
masers (Tofani et al 1995, Minier et al 2001, Edris et al 2005) and an outflow/jet system 
traced in CO, SiO, H 2 and cm- wavelength continuum emission. It has a rotating disk of 
radius ~ 5000 AU, traced in ammonia and CS (Zhang, Hunter & Sridharan 1998; Cesaroni 
et al 2005) and a possible smaller ~ 1000 AU radius disk inferred from interferometric 
CH 3 CN measurements (Cesaroni et al 1999). Sensitive cm-wavelength continuum images 
have uncovered emission south of the central object suggesting the presence of a companion 
(Hofner et al 1999), also indicated by the precession of the H 2 jet from this source (Shepherd 
et al 2000, Cesaroni et al 2005). We note that IRAS20126+4104 is a member of a large 
sample of HMPO candidates we have been studying systematically (Sridharan et al 2002, 
Beuther et al 2002a, b, Williams, Fuller & Sridharan 2004, 2005). 

Here, we report the infrared detection of a compact disk/binary system in IRAS20126+4104, 
consisting of a small scale circumstellar disk and a secondary source. 
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2. Observations 

The observations were carried out during 15-18 Aug, 2000. We imaged IRAS20126+4104 
using the near-infrared cameras UFTI (1-2.5 /im) and TUFTI/IRCAM (1-5 /im) on the 
United Kingdom Infrared Telescope (UKTRT). 1 TUFTI is a 256x256 InSb Array with a 
field of view of 20.7" and 0.08" pixels. UFTI has a 1024x1024 HgCdTe Array with 0.09" 
pixels and a field of view of 93". The seeing during these measurements was excellent, at 
~ 0.3" (K band). Observations were made with broadband K, L', and M' filters (central 
wavelengths of 2.2, 3.8 & 4.7/im). The K band observations employed a 11" jitter and the 
L' and M' observations used 8.5", 7.5" and 6.5" jitters and a 2" jitter with a 30" chop. 
Standard data reduction steps were used to obtain the final images. Several calibration stars 
were observed for reliable photometry. 

We used K-band objects from the 2MASS survey to fit separate plate solutions for the 
K and L' images. The formal errors for the plate solutions are 0.15" rms, consistent with the 
2MASS astrometric accuracy 2 . The error on the relative registration is 0.1" rms, measured 
using the positions of the same objects in the independently registered K and L' images. 
Since the M' image does not include any reference objects, we cross-correlated it with the L' 
image to find the best position offset. We estimate the absolute registration uncertainty to 
be ~ 0.15". 



3. Results and Discussion 

As shown in Figure 1, we detect two, previously unknown, faint lobes of emission in the 
K-band, K-SE & K-NW, near the center of the IRAS20126+4104 outflow region. Previous 
images only showed the bi-polar features on larger ~ 10" scales oriented at about the same 
position angle as the molecular and the cm-wavelength outflows (Cesaroni et al 1997, Ayala 
et al 1998, Shepherd et al, 2000). The new emission objects are separated by a dark lane 
with the line joining them oriented NW-SE, similar to the larger scale bi-polar features. 
Figure 2 (left panel) shows the L' and K emission in the central 2", where the center of the 
dark lane is seen to nearly coincide with the location of the mm-wavelegth emission from 
interferometric images. The L' emission over the region lies close to the dark lane seen in 
the K band. The peak of this central emission is nearly centered on the dark lane, within 



1 Thc UKIRT is operated by the Joint Astronomy Centre on behalf of the U.K. Particle Physics and 
Astronomy Research Council 

2 2MASS Explanatory Supplement, http://ipac.caltech.edu/2mass/releases/allsky/doc/expl-sup.html 
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registration errors. Emission is also seen towards the K-band lobes. In addition, there is a 
point source like object seen ~ 0.5" southwest of the center. This second object is at one 
end of, and in line with, the dark lane. The center panel in Figure 2 shows the M' emission 
contours superposed on the K image. The M' emission has a similar morphology as that in 
the 17 band, but with the central object more prominent. The correspondence between the 
emission in the L' and the M' bands is shown in the right panel. The source detected at L' 
to the southwest of the peak is also seen in the M' image, confirming its reality. This also 
validates our method of cross-correlating L' and M' band images for registration, which is 
sensitive only to the strongest emission features in the images. 

Bi-polar infrared morphologies similar to the ones found in our data are frequently seen 
in near infrared observations of low mass young stars, for example the HST images of T Tauri 
stars by Padgett et al. (1999). Towards the low mass objects these morphologies have been 
interpreted as due to the presence of nearly edge-on circumstellar disks. Therefore we suggest 
that our images represent the first detection of a small scale disk towards a young high-mass 
star, as a silhoutte at NIR wavelengths. The width of the dark lane, which corresponds to 
the thickness of the disk is ~ 0.5" and the diameter is less certain at ~ 1". The implied linear 
dimensions of the disk are ~ 850 AU thickness, for radii < 1000 AU. This is comparable to 
the size of the disk inferred from the CH 3 CN data of Cesaroni et al (1999). These sizes are 
significantly larger than the small-scale disks seen towards low-mass star-forming regions as 
may be expected for higher-mass proto-stars. 

The presence of the second source makes this object more interesting. As seen in Figs 
2-4, this second object is located on the same side of the central object as the cm- wavelength 
jet of Hofner et al (1999; shown by a cross mark to the south). Given the closeness of the 
two L' & M' band objects, (~ 0.5" = 850 AU), we may be seeing a young high-mass binary 
system. We suggest that the second object to the southwest is the companion responsible 
for the precessing H 2 emission outflow (Shepherd et al 2000, Cesaroni et al 2005). It is 
noteworthy that the companion is in line with the disk. It is probably close to the equatorial 
plane, although due to unresolved projection effects, this cannot be said with certainty. We 
point out that in the CH 3 CN disk inferred by Cesaroni et al (1999), the emission peaks 
appear to cut off beyond the location of the companion. The location of the companion does 
not match the location of the cm-jet, which raises the possibility of a triple system. 

3.1. Extinction Estimates 

We now make some very simple and approximate estimates of the extinction due to 
the envelope surrounding the disk-binary system and that due to the disk, separately. We 
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use the observed colours of the binary and the bi-polar features along with the interstellar 
reddening law from Koorneef (1983) (in particular, E L _ M /A V = 0.019 & E K _ L /A V = 0.045) 
for this purpose. For all stellar objects, intrinsic K-L' and L'-M' colours of ~ are assumed. 
Two cases, scattered and direct, are considered, where line of sight reddening occurs with 
and without preceding Rayleigh scattering (A~ 4 dependence; K-L' = —2.4 & L'-M' = —1.0). 

The magnitudes and colours, from aperture photometry, for the whole region, the central 
lobes and the binary components (central and companion) are presented in Table 1. For the 
K band, using the fluxes of the H 2 knots (Cesaroni et al 1997) and spectra from Ayala et 
al (1998), a line emission contamination of a few percent or < 0.1 magnitude is estimated. 
Contamination in the L' and M' bands are not considered. This is reasonable for the L' band 
due to its large bandwidth. The integrated emission over the whole region, excluding obvious 
point-like sources, is in agreement with the K band flux reported by Cesaroni et al (1997). 
The colour of this emission is consistent with light from main sequence stars, implying 120 & 
116 magnitudes of Ay (mean 118) for the scattered case, and ~ 65 magnitude for the direct 
case. Therefore, an Ay of 65 - 118 is considered to be the extinction due to the envelope. 

For the central lobes, the extinctions (scattered case) implied are 113 and 127 (mean 
120) from K-L' colours, similar to the upper limit estimate for the ~ 10" features. Since the 
emission from central and companion objects and the bi-polar lobes cannot be well seperated, 
the L' magnitudes are lower limits (the lobes may be fainter), and the extinctions, upper 
limits. In summary, the envelope extinction is in the range 65 - 120. 

Turning to the binary itself, assuming ~ for the intrinsic colours of the stars, the 
limits on the observed K-L' colours imply a lower limit (not detected in the K band) to Ay 
towards the central and the companion objects of ~ 96 & 91, for the direct case and 149 & 
144 for the scattered case. From the L'— M' colours, the values are 168 & 148 and 221 & 200, 
for the two cases, consistent with the lower limits from K-L' colours. The ranges 168-221 
and 148-200 will be used for the two objects for further discussion. This translates to a H2 
column density of ~ 2xl0 23 cm" 2 (N H JA V = 0.94 x 10 21 ; Frerking, Langer & Wilson, 1982). 
In comparsion, Beuther et al (2002a) estimated an extinction of 553 magnitudes for the dust 
core (11" beam), or 276 magnitude up to the center. 

It is interesting that the extinction towards the companion is less than that towards the 
central object: we speculate that this may be due to a gap cleared by the companion outside 
of which the material may be less dense. As noted before, the CH 3 CN emission peaks in 
Cesaroni et al (1999), which trace high densities, are distributed within the area bounded 
by the companion which supports such a picture. 
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3.2. Magnitudes and Masses 

These measurements allow an estimate of the range of masses for the disk and NIR 
magnitudes for the embedded stellar objects. With an envelope extinction of ~ 65-120, the 
excess extinctions towards the two objects are in the range A v ~ 48-166 and ~ 28-135. 
This extinction may be attributed to the disk and translates to N H2 ~ 10 23 cm -2 . For 
a 1000 AU x 2000 AU disk, a mass of ~ 0.1 M is indicated. This is quite uncertain 
due to the many assumptions, and our inablility to reliably separate the contributions to the 
measured magnitudes from the two embedded young stars and the bi-polar lobes. Additional 
uncertainties are due to part of the L' and M' band emission coming from hot dust in the disk 
which would change the colours, and the mass of such a component is also not considered. 
Flaring of the disk and its inclination to the line of sight, which would increase the mass 
estimates, are not considered. Therefore, we beleive 0.1 M to be an approximate lower 
limit to the mass within 1000 AU from the central star. In comparsion, Cesaroni et al (2005) 
estimate the mass of the disk within 5000 AU from the star to be 1-4 M . They estimated 
peak H 2 column densities for the disk of 1.5-3.4xl0 22 cm -2 from C 34 S measurements (using 
numbers listed in their Table 4) compared with our ~10 23 cm~ 2 . We note that our images do 
not constrain the radius of the disk - it can be larger than ~1000 AU, in which case our mass 
estimate would go up. Also, observations such as ours selectively probe the lowest extinction 
paths through which photons escape, leading again to lower estimates for extinctions and 
masses. We note that Fuller, Zijlstra & Williams (2001) estimated a similar low mass for a 
proposed NIR disk in the high-mass jet object IRAS18556+0136. 

What are the masses implied for the two young stars? Assuming main sequence objects 
and total extinctions of Ay ~ 168-221 and ~ 148-200 for the two, the extinction corrected 
L' magnitudes are: 3.6 - 1.3 and 4.7 - 2.4 (A L /A V = 0.045; from Koorneef, 1983). With a 
distance modulus of 11.2, the absolute L' magnitudes implied are: —7.6 — 9.9 & —6.5 — 8.8 
corresponding to very early O spectral types - 04 and earlier (Wegner 2000; Serabyn, Shupe 
& Figer, 1998). For the primary object, Cesaroni et al (2005) derived a mass of ~ 7 M 
using the Keplerian velocity field they observed, singnificantly changing earlier estimates of 
20-24 M and implied ZAMS spectral types of 08-09 (Cesaroni et al 1999, Zhang, Hunter & 
Sridharan, 1998). Again, we emphasize that our estimates are to be considered approximate 
- in addition to the previously mentioned uncertainties, any emission in the L' and M' bands 
from the disk would lower the magnitudes estimated for the stars - and it is almost certain 
that the stars are not as early types as implied. Also, the main sequence assumption is 
questionable and the absolute NIR magnitudes of the brightest main sequence stars remain 
poorly characterized. The companion is probably only ~ 1 magnitude fainter, which makes it 
significantly less massive, but comparable. A significant mass for the companion will appear 
to be in conflict with the recent observations of Keplerian velocity profile (Cesaroni et al 
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2005, Edris et al 2005). However, it is not clear if the Keplerian profile holds good close to 
and inside the 0.5" separation of the companion. Higher resolution observations are critical 
to resolving the disagreements. Nevertheless, our numbers do indicate that we are indeed 
dealing with a high-mass system. 

Given the marginal linear resolution of the images, we have refrained from detailed 
modelling which should wait until higher resolution images are obtained with adaptive optics. 
Fortunately, there is a bright star ~ 12" from the center which makes IRAS20126+4104 an 
ideal candidate for such studies. 



4. Conclusions 

We presented high resolution NIR images of IRAS20126+4104. The K-band image 
shows a dark lane separating two, new, l"-scale bi-polar emission features, with the emission 
peaks in the longer wavelength L', M' and mm-wave bands lying close to the dark lane. This 
is interpreted as due to a nearly edge-on disk of 850 AU thickness for radii < 1000 AU, imaged 
and measured for the first time at infrared wavelengths. A new, point-like, second object, 
found to the southwest of the central object, and in line with the disk, is interpreted as a 
binary companion with a separation of ~ 1000 AU. Estimates of the NIR magnitudes of the 
objects, although very approximate, indicate a high-mass system. Through their striking 
resemblence to similar low-mass systems, the new images demonstrate that the processes 
involved in the formation of low and high mass stars may have much in common: viz., 
multiple star formation in an accretion disk setting. Further higher resolution studies will 
be extremely rewarding in modelling the system and therefore to address a number of issues 
related to high-mass star-formation. 

We thank the referee, Riccardo Cesaroni, for careful comments which significantly im- 
proved the paper. 
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Fig. 1. — K-band image of IRAS20126+4104 (J2000 co-ordinates). The central region, shown as 
a box, is enlarged in the right panel, where two faint emission lobes are seen. 
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Fig. 2. — K, L' and M' emission from the central region, left: K-band image L'-band emission 
contours. The asterisks mark water masers (Tofani et al 1995). The central cross mark locates the 
mm- wave emission (an average of the peak positions at 1.3, 3 and 7-mm), also coincident with one 
peak of the VLA 3.6-cm emission. The cross mark to the south is the location of the second 3.6cm 
emission jet (Hofner et al 1999). center: K-band image with L'-band contours, right: L'-band 
image with M'-band contours. 
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Table 1. Table 1 - NIR Magnitudes, Colours and Extinctions 



Object 


Magnitude 




Colour 


Ay t direct 


^4d, scattered 




K 


L' 


M' 


K-L' L' 


-M' 


(K-L') 


(L'-M') 


(K-L') 


(L'-M') 


Whole 


7.7 


4.7 


3.5 


3.0 


1.2 


67 


63 


120 


116 


SE lobe 


13.5 


10.8 




2.7 




60 




113 




NW lobe 


14.9 


11.6 




3.3 




73 




127 




Central 


>15.5 


11.2 


8.0 


>4.3 


3.2 


>96 


168 


>149 


221 


Companion 


>15.5 


11.4 


8.6 


>4.1 


2.8 


>91 


148 


>144 


200 



Note. — The direct coulmns assumed extinguished star light from stars with intrinsic 
K-L' and L'-M' colours of 0. The scattered columns assumed Rayleigh scattering before 
extinction. 
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